Quantum computing using dissipation 



A Beige 

Blackett Laboratory, Imperial College London, Prince Consort Road, London, SW7 2BW, 

m ■ UK 

o 
o 

' Abstract. The principal obstacle to quantum information processing with many qubits is 

decoherence. One source of decoherence is spontaneous emission which causes loss of energy 
and information. Inability to control system parameters with high precision is another possible 
source of error. As a solution we propose quantum computing experiments using dissipation 
based on an environment-induced quantum Zeno effect. As an example we present a simple 
scheme for quantum gate implementations with cold trapped ions in the presence of cooling. 
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^ ■ 1. Introduction 

\Q \ Following the theoretical formulation of quantum computing [ 1 1 and the first algorithms 
for problems which can be solved more easily on a quantum computer than on a classical 
computer 10 13 the practical implementation of" such a device has become a challenging task. 
Initial steps have already been taken. Quantum bits (qubits) can be realised for instance by 
storing the information in a superposition of the internal states of two-level atoms. However 
building systems with many coupled qubits remains a huge challenge. Many demands must be 
met: reliable qubit storage, preparation and measurement, gate operations with high fidelity 
and low failure rate and scalability of the system to many qubits. The biggest problem is 
q-i posed by decoherence which can cause the loss of information. 

In this paper we review the idea of quantum computing using dissipation |@J|3 which 
might help to overcome the decoherence problem in some systems. On the contrary, it might 
even be useful to introduce an additional spontaneous decay channel into a system in order to 
create a realm of new possibilities to implement gate operations between qubits. To illustrate 
this we present a scheme for the realisation of quantum computing between cold trapped ions 
in the presence of cooling. The ions are stored inside a linear trap and each qubit is obtained 
from two different ground states of the same ion, as in |6||7|. 

In the quantum computing scheme proposed here, the role of cooling during gate 
operations is twofold. On one hand it decreases the sensitivity of the scheme with respect to 
heating [ 8 1 . On the other hand, the presence of the cooling lasers introduces a decay channel 
into the system whose presence leads to a restriction of the time evolution of the system onto 
the computational subspace and facilitates the possibility for quantum gate implementations 
within one step. The probability for photon emission is relatively small since the system 
remains to a good approximation within a decoherence-free subspace [0 [TUJ [TT|. If heating 
nevertheless populates the vibrational mode or gate failure moves the system out of the 
decoherence-free subspace, then photons are emitted at a high rate JT2l 151 IT3l . This can 
be detected and the computation can be restarted. 

Many schemes for quantum computing with trapped ions have already been proposed. 
Some of them require cooling of the ions into the ground state of a common vibrational 
mode. That this is possible has been demonstrated recently in Innsbruck iflUl . where the 
Cirac-Zoller controlled-NOT quantum gate [15| has been implemented with the help of six 
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concatenated laser pulses individually addressing each of the two ions. At the same time, 
the group in Boulder demonstrated a robust, high-fidelity geometric two ion-qubit phase gate 
in the laboratory. This was achieved with a sequence of laser pulses and without individual 
addressing of the ions ifTol . But are these schemes really suitable for quantum computation 
with many qubits? Finding reliable ways to scale present schemes to many qubits might 
require further simplifications of the experimental setup without decreasing the precision of 
gate operations. 

2. Non-Hermitian Hamiltonians and no-photon time evolutions 

Before we discuss the coherent control of an open quantum system with the help of dissipation 
in more detail, let us first introduce the theoretical model for describing the time evolution of 
the system under the condition of no photon emission. We give a short review of the quantum 
jump method 1T7l which is equivalent to the Monte Carlo wave-function [18] and the quantum 
trajectory lTl9l approaches ll20l . 

2.1. The quantum jump approach 

In the last decades, several quantum optics experiments have been performed studying the 
statistics of photons emitted by single quantum mechanical systems, like one or two laser- 
driven trapped atoms or ions. Effects have been found that would be averaged out in the 
statistics of photons emitted by a whole ensemble and which cannot be predicted with the help 
of expectation values calculated for a statistical ensemble lETl . New formalisms describing 
single systems interacting with the environment had to be developed. 

A typical example for such an experiment is electron shelving [22|, i.e. the occurance of 
stochastic macroscopic light and dark periods in the resonance fluorescence of a laser-driven 
atom with a metastable state. Another example is the two-atom double-slit experiment by 
Eichmann et al. ll23l which demonstrated that the photons emitted by two atoms at a fixed 
distance can create an interference pattern on a distant screen. These experiments suggest that 
the effect of the environment on the state of the atoms is the same as the effect of rapidly 
repeated measurements of whether a photon has been emitted or not ll24l 1251 . From this 
assumption the quantum jump approach has been derived. 

Assume that a measurement is performed on a quantum optical system surrounded by 
a free radiation field initially in its vaccum state |0 p h) and prepared in \ip) determining after 
a time At whether or not a photon has been created. If H is the Hamiltonian including the 
interaction of the system with its environment, the state of the system equals 



under the condition that the free radiation field is still in the vacuum state. For quantum 
optical experiments, it has been shown that the dynamics under the conditional time 
evolution operator U cont &(At,0), defined by the right hand side of ([T]), can be summarised 
by a Hamiltonian H cond that is largely independent of the choice of At. The conditional 
Hamiltonian H cond is non-Hermitian and the norm of a state vector developing with H cond 
decreases in general in time. Under the condition of no photon emission, the state of the 
system equals at time t 



ph ) C/cond(At,0)|^) = |0 ph )(0 ph | U(At,0) |0 ph )|V> 



(1) 




(2) 



(3) 
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Figure 1. Level scheme of the two ions involved in the gate operation. Each qubit is obtained 
from the ground states |0) and |1) of one ion. In addition a metastable state |2), a rapidly 
decaying state |3) with decay rate and two strong laser fields with coupling strength 
gj = | rjjCtj and detuning v are required. 



The non-Hermiticity of the conditional Hamiltonian i/ CO nd and the continuous decrease 
of the amplitude of the unnormalised state vector U con d(t, 0) \ip) reflect that the observation of 
no photons reveals information about the system. The longer no photon is emitted the more 
unlikely it becomes that there is excitation that might cause an emission and the amplitudes 
of states with spontaneous decay rates decrease exponentially [|2"6Tl. 

2.2. Example: Two four-level ions inside a linear ion trap 

The level scheme of the two cold trapped ions that we consider in this paper in order to 
discuss the role of dissipation in quantum gate realisations is shown in Figure [T] Each qubit 
is obtained from two different atomic ground states |0) and 1 1) of the same ion. In addition, 
a metastable state |2) and a rapidly decaying level 3 are required. First, the ions have to be 
cooled into the ground state of a common vibrational mode. Two strong laser fields detuned by 
the frequency v of a common vibrational mode should be applied. The laser field coupling to 
the 1-2 transition establishes a coupling between the two qubits involved in the gate operation. 
The laser field driving the 1-3 transition represents the laser cooling setup and can be replaced 
by any other laser cooling configuration without changing the effective time evolution of the 
system. 

In the following, we denote the spontaneous decay rate of level 3 by T 3 while b and 
are the annihilation and creation operator of a phonon in the common vibrational mode. The 
coupling constant of this mode to the atomic 1-j transition equals gj = | TfjQj where Qj is 
the Rabi frequency of the applied laser field and rjj is the Lamb Dicke parameter depending 
on the characteristics of the ion trap. Proceeding as in ifTTll . one finds that the conditional 
Hamiltonian within the dipole and the rotating wave approximation and in the interaction 
picture with respect to the free Hamiltonian equals 

2 2 

#cond = X>[<?2 |l>i(2| & f + 93 & f - h.c. ] - £ i%T 3 |3) l (3| . (4) 
i=i i=i 

Here the Lamb-Dicke regime and the condition v ^> Q2, has been assumed, as in lfT31 . 
3. Quantum computing using dissipation 

The basic idea of quantum computing using dissipation is to utilise the coherent no-photon 
time evolution given by H conA f° r the implementation of gate operations. Whenever a photon 
emission occurs the computation fails and has to be repeated. Nevertheless, we show that 
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fidelities F and success rates P close to one can be achieved for a wide range of experimental 
parameters. The reason that spontaneous emission from the ions is in the following negligible 
is that the system remains during the whole computation in a decoherence-free (DF) state 
El[T0l[TT]|. In this section we define DF states and describe a mechanism which restricts the 
time evolution of a system onto its decoherence-free subspace (DFS). 

3.1. Decoherence-free states 

The DFS of a system is a subspace of states whose population does not lead to decoherence. 
Using the quantum jump approach, a state is DF if Po(i, ip) — 1 for all times t ifTTI . Hence, 
the DFS is spanned by the eigenvectors of the conditional Hamiltonian with real eigenvalues 
Aj. The eigenvectors |Aj) of H con d are in general non-orthogonal. It is therefore useful to 
introduce the reciprocal basis vectors |A J ) with (A J |Aj) = 5ij and to write the conditional 
Hamiltonian as 

-ffcond — /]\ • (5) 

i 

Suppose that all non-DF states couple strongly to the environment and populating them leads 
typically to a photon emission within a time At. The time At is greater than a certain minimal 
size which can be determined from the quantum jump approach. Provided that the eigenvalues 
Afc corresponding to non-DF states fulfil the condition 

e -iA fc At/ft = 0) (6) 

the no-photon time evolution operator becomes ll271 

*7 CO nd(A*,0) = Yl e- iX > At / h |A«> <A*| . (7) 

i:|Ai>eDFS 

This operator projects every initial state onto the DFS. The action of the environment over a 
time At can therefore be interpreted as a measurement whether the system is DF or not. The 
probability for no emission in At equals the probability to be in a DF state. 

3.2. An environment-induced quantum Zeno effect 

These continuous measurements caused by the environment lead to a realm of possibilities to 
manipulate a system within the DFS. Actually, any arbitrary interaction can be used as long as 
its typical time scale is much longer than At defined by condition © [4|. This restriction of 
the system onto the DFS can intuitively be understood with the help of the quantum Zeno 
effect ll28l . Within At, a weak interaction can only transfer population proportional At 
out off the DFS. During the next measurement the system is found in a non-DF state with 
a probability proportional At 2 . Otherwise it is projected back onto the DFS by the time 
evolution operator Q). Important is that the probability to find the system always in a DF 
state, i.e. T/At times if T is the gate operation time, goes in the limit of weak interactions, 
i.e. for At/T -> 0, to one. 

However, the time evolution of the system inside the DFS is not inhibited. Using first 
order perturbation theory with respect to the weak interaction, © and the assumption that the 
system is initially in a DF state, one can show that the conditional time evolution in At is the 
same as the one that one obtains from the Hamiltonian 

H e tt = IP DFS Hcond -PdFS (8) 

using the same approximations. Here IPbfs = X^ia^gdfs |<^t)(^*| ls me projector onto the 
DFS. The effective Hamiltonian H c s is used in the following to find the appropriate laser 
configuration for the realisation of certain gate operations. 
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One motivation among others to use an environment-induced quantum Zeno effect for 
quantum computing is the simplicity of the resulting schemes. As we see in the next section, 
the DFS of the cold-ion system contains in addition to all ground states, highly entangled 
states. Through populating these states entanglement between qubits can be created during 
the effective time evolution even if this is in general not possible using only a single laser 
pulse. 

3.3. Scaling of probabilistic quantum computing schemes 

Let us now estimate the probability of finding the result of a whole computation assuming that 
each gate can be performed with maximum fidelity (as it applies as long as the corrections 
to the desired effective time evolution are not too big 0) but only with a finite success rate 
Pq. The probability of implementing an algorithm of iV gates faultlessly is P$ and grows 
exponentially with N. On the other hand, if one always knows whether an algorithm has 
failed or not, the computation can be repeated until a result is obtained. The probability for 
not having a result after M runs equals 

-Pno result = ( 1 — Pq ) ■ (9) 

For large iV this is approximately exp ( — M P^). Many repetitions might be necessary to 
implement a computation. However, for smaller numbers of iV and if P is sufficiently close 
to unity, the failure probability is already nearly negligible for M ~ iV '. For example, if 
Po = 95 % and an algorithm with N = 50 gates is performed, then repeating the computation 
50 times yields a success rate above 98 %. 

4. Quantum computing with cold trapped ions in the presence of cooling 

We now consider again the concrete setup described in subsection 12.21 and describe the 
realisation of two-qubit quantum gates between cold trapped ions in the presence of cooling 
with the help of additional weak laser fields. 

4.1. Decoherence-free states of the ions 

To predict the effect of these fields we first determine the decoherence-free (DF) states of the 
system in the absence of any additional interaction. As shown in subsection 13.11 they can 
be calculated by finding the eigenvectors of the Hamiltonian © with real eigenvalues. This 
yields that the decoherence-free subspace (DFS) of the system contains only superpositions of 
states with the ions either in the state 1 00) combined with an arbitrary state of the vibrational 
mode or both ions in 1 1 ) , 1 1 0) , 1 1 1 ) or in the antisymmetric and maximally entangled state 

|o> = ^[|12>-|21>] (10) 

while the vibrational mode is not populated. Here all real eigenvalues are zero and the system 
does not evolve as long as no additional interaction is applied. 

4.2. Single laser pulse quantum gates 

To realise gate operations between the two qubits formed by the ground states of the ions 
weak laser fields are applied in addition to the strong lasers shown in Figure [T] Let us denote 

the Rabi frequency of the laser with respect to the j-2 transition in ion % by ilj and assume 
that 

nS°<<fe, 93 and T 3 . (11) 
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Figure 2. Success rate and fidelity under the condition of no photon emission of a single 
CNOT gate as a function of the Rabi frequency fl for T3 — 2^/37 g2, 33 = V2 52 and for the 
initial qubit states |00) (a), 1 10) (b), |11) (c) and [ 1 10) - |11) ]/y/2 (d). The gate success rate 
and fidelity are always maximal if the ions are initially prepared in |01). 



The conditional Hamiltonian H cond of the system then becomes the sum of the Hamiltonian 
© and the laser Hamiltonian 

#lase r = E E imf\j)i(2\+h.C. (12) 
i=l,2 j=0,l 

If the time evolution of the system is restricted onto the DFS, as predicted in the previous 
section, then it is effectively given by the Hamiltonian ® which equals 

H eS = ^h{- |01) + Q { o ] 1 10) + (SI? - SI?) |11) ] (a\ + h.c. (13) 

This Hamiltonian can be used to implement quantum gate operations if the operation time 
T is chosen such that at the end of each gate again only qubit states are populated. A more 
detailed analysis of the no-photon time evolution of the ions is given in 0. 



4.3. The controlled-NOT gate 

As a concrete example we consider the controlled NOT (CNOT) gate. If ion 1 contains the 
target qubit and ion 2 provides the control qubit this gate corresponds to the time evolution 
operator 

^gate = |00)(00| + |01>(01| + |10>(11| + |11)(10| . (14) 
The easiest way to realise a CNOT gate is to couple one laser with the (real) Rabi frequency 
SI = fij to the 1-2 transition of ion 1 and another one with the same Rabi frequency to the 
0-2 transition of ion 2 which yields the effective Hamiltonian 

#eff=2^ftft[|10)-|ll)](a|+h.c. (15) 

If the duration of the laser pulse equals T = 2n /SI the resulting time evolution can be shown 
to be exactly the desired operation 0. 

Figure El results from a numerical solution of the no-photon time evolution given by the 
sum of the Hamiltonians © and (fT2l . For very small Rabi frequencies, the gate success 
rate and fidelity is for all initial states close to one. For larger values of f2, the P (T,ip) 
decreases and is for = 0.2 g 2 as low as 73%. The gate fidelity is in this case still above 
98.4 %. The smallest gate success rate is found when the atoms are initially in |10), |11) or in 
a superposition of these two states. In general, success rates P > 90 % are achieved as long 
as n< 0.07 # 2 - 
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Figure 3. Fidelity for a single CNOT gate under the condition of no photon emission as a 
function of the spontaneous decay rate T^/ 92 for 53 = T^, f2 = 0.01 32 and for the initial 
qubit states |00) (a), 1 10) (b) and [ 1 10) - 1 11) ]/V2 (c). If the ions are initially prepared in 
1 01), the gate success rate equals one. For \ip) — 1 10) the fidelity is about the same as in graph 
(c). 



Finally, we would like to comment on the role of the spontaneous decay rate of level 
3 in the scheme. Figure |3] shows the fidelity of a single CNOT gate as a function of T 3 
and for g 3 = T 3 . In the chosen parameter regime, the effective damping rate of unwanted 
population in non-DF states can be shown to equal gl/T 3 = T 3 to a very good approximation 
and the effective decay rate of non-DF states increases linearly in T 3 . Figure |3] confirms that 
the presence of the auxiliary dissipation channels is crucial for the scheme to work! For 
small damping rates the mechanism which restricts the time evolution of the system onto the 
computational subspace fails and the minimum gate fidelity is well below 50 %. 

5. Conclusions 

The paper reviews the idea of quantum computing using dissipation which offers a great va- 
riety of possibilities to implement gate operations in the presence of only one dissipation 
channel in a system. As an example we described the possibility to implement precise quan- 
tum gates between cold trapped ions in the presence of cooling of a common vibrational mode 
and within one step. Because of its simplicity and robustness against parameter fluctuations 
[7 1, the proposed scheme might help to increase the number of qubits in present quantum 
computing experiments. 
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